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ABSTRACT 

A new model of the observed universe, using solutions to the full Einstein equations, is developed from 
the hypothesis that our observable universe is an underdense bubble, with an internally inhomogeneous 
fractal bubble distribution of bound matter systems, in a spatially flat bulk universe. It is argued 
on the basis of primordial inflation and resulting structure formation, that the clocks of the isotropic 
observers in average galaxies coincide with clocks defined by the true surfaces of matter homogeneity 
of the bulk universe, rather than the comoving clocks at average spatial positions in the underdense 
bubble geometry, which are in voids. This understanding requires a systematic reanalysis of all observed 
quantities in cosmology. I begin such a reanalysis by giving a model of the average geometry of the 
universe, which depends on two measured parameters: the present matter density parameter, f2 m , and 
the Hubble constant, H Q . The observable universe is not accelerating. Nonetheless, inferred luminosity 
distances are larger than naively expected, in accord with the evidence of distant type la supernovae. 
The predicted age of the universe is 15.3 ± 0.7 Gyr. The expansion age is larger than in competing 
models, and may account for observed structure formation at large redshifts. 

Subject headings: Cosmology: theory - Cosmology: large-scale structure of universe — Cosmological 
parameters — Cosmology: early universe [ arXiv: gr-qc/0503099 ] 


1. INTRODUCTION 

Observations in the past decade have been interpreted 
as suggesting that 70% of the matter-energy density in the 
universe at the present epoch is in the form of a smooth 
vacuum energy, or “dark energy”, which does not clump 
gravitationally. This is supported by two powerful inde¬ 
pendent lines of observation. Firstly, type la supernovae in 
distant galaxies (Perlmutter et al. 1998, 1999; Riess et al. 
1998, 2004) are dimmer than would be expected in stan¬ 
dard Friedmann-Robertson-Walker (FRW) models, espe¬ 
cially when it is noted that many independent dynamical 
estimates of the present clumped mass fraction, U m , sug¬ 
gest values of order 20-30%. Secondly, observations of the 
power spectrum of primordial anisotropies in the cosmic 
microwave background radiation (CMBR), most recently 
by the WMAP satellite (Bennett et al 2003), indicate that 
the universe appears to be spatially flat on the largest of 
scales. 

A cosmological constant, or alternatively dynamical 
dark energy, is most commonly invoked to explain the ob¬ 
served cosmological parameters, even though a fundamen¬ 
tal origin for such dark energy remains one of the pro- 
foundest mysteries of modern physics. However, even in 
the presence of dark energy at the present epoch, a number 
of problems remain. One of the most significant problems 
is that the epoch of reionization measured by WMAP ap¬ 
pears at a redshift of order z~20lg , indicating that the 
first stars formed much earlier than conventional models 
of structure formation would suggest. The detection of 
complex galaxies at relatively large redshifts compounds 
the conundrum (see, e.g., Cimatti et al. 2004, Glazebrook 
et al. 2004). 

In recent work, Kolb, Matarrese, Notari and Riotto 
(2005) have proposed a profoundly different resolution of 
the “cosmological constant problem”. They reason that 


primordial inflation will have produced density perturba¬ 
tions many times larger than the present horizon volume. 
We should not view the observable universe as typical of 
the universe on scales larger than our particle horizon, and 
the values of cosmological parameters observed should not 
necessarily be taken as typical of the whole. Furthermore, 
they suggest that our present observations might be com¬ 
patible with the observed universe being an underdense 
bubble in an otherwise spatially flat k = 0 FRW universe 
with an energy density fl TOT = 1 in ordinary matter. 
Such assumptions are indeed supported by detailed numer¬ 
ical and analytic calculations made in inflationary models 
a decade ago (Linde, Linde and Mezhlumian 1996). 

The possibility of inhomogeneously defined clock rates is 
commonly accepted in studies of inhomogeneous cosmolo¬ 
gies, such as the Lemaitre-Tolman-Bondi (LTB) models. 
(For a review see Krasinski (1997).) In this Letter , I will 
argue, however, that primordial inflation can give rise to a 
particular structure in inhomogeneous cosmologies, which 
allows for a homogeneous cosmic time on the scales of the 
first bound systems which form, which differs from the “co¬ 
moving” time parameter of the average late epoch geom¬ 
etry. This novel feature is what distinguishes the present 
model from previous studies of inhomogeneous cosmolo¬ 
gies, as it leads to a new solution of the fitting problem 
(Ellis and Stoeger 1987). 

2. OBSERVERS, CLOCKS AND THE FITTING PROBLEM IN 
INFLATIONARY COSMOLOGY 

It is a consequence of the inflationary paradigm that 
a spectrum of initially small density perturbations is 
stretched to all observable scales within our past light cone, 
and also to scales beyond our particle horizon. The fact 
that this is true for the past light cone is well supported 
by the CMBR. The hypothesis that such perturbations 
should extend to super-horizon scales is a feature of most 
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inflationary models, independent of their details. 

One important realisation is the fact that since primor¬ 
dial inflation ended at a finite very early time, the scale 
of super-horizon sized modes, although huge, must have a 
cut-off at an upper bound. Beyond that scale we assume 
the universe is described by a spatially flat bulk metric 

dsLik = ~dr 2 + a 2 (r) (dx 2 + d y 2 + dz 2 ), ( 1 ) 

where a(r) = ai(r/T i ) 2 /[ 3 ( 1 +t< ’)l, and we use units with 
c = 1. The precise bulk equation of state Pbuik = rupbuik 
proves to be inconsequential if w > — 1 , though in accord 
with the principles of the model we assume only ordinary 
matter and radiation. 

Although the observable universe will undoubtedly be 
embedded in many regions of under- and over-density, like 
the smallest figure inside a Russian doll, it is nonethe¬ 
less reasonable to assume that provided the density per¬ 
turbation immediately containing our observed universe 
extends sufficiently beyond our horizon then there is a 
super-horizon sized underdense bubble containing the ob¬ 
servable universe, with matter density equal to the aver¬ 
age matter density we measure, which we can model as a 
super-horizon sized underdense region, S. 

The matter distribution inside S is assumed to be inho¬ 
mogeneous and ‘‘fractal’’, in accord with calculations from 
primordial inflation (Linde, Linde and Mezhlumian 1996) 
and in accord with all the evidence of observations of the 
actual universe, with its large-scale structure of stream¬ 
ing motions of clusters of galaxies, bubbles and voids. 
Nonetheless, despite this inhomogeneity we do observe an 
average isotropic Hubble flow. Thus at some level the first 
step in solving the fitting problem (Ellis and Stoeger 1987) 
must involve the approximation of the inhomogeneous ge¬ 
ometry of S by an average spacetime which depends on 
a single cosmic scale factor, viz., the spatially open FRW 
geometry 


dS 2 = -At 2 


a 2 (t) 


dr 2 


1 + 7"' 


+ r 2 (dd 2 + sin 2 6 A(j) 2 ) 


( 2 ) 


where parametrically in terms of conformal time, 77 , 


a = -- —(cosh 77 — 1 ), 


Hit = 


2(1 - no 

2(1 — fb) 3/2 


(sinh 77 — 77 ). 


( 3 ) 


where f \ is an initial density parameter of the bubble, S , 
and a* a constant to be set. 

Usually the first step after writing down the metric (2) is 
to identify comoving observers in this geometry with ide¬ 
alised isotropic observers, namely observers who measure 
no dipole anisotropy in the CMBR. One identifies typical 
stars in typical galaxies with such observers, apart from 
the small effects of peculiar velocities. However, even in 
an homogeneous isotropic FRW model this implicitly as¬ 
sumes a solution of the fitting problem in terms of a com¬ 
plicated matching of asymptotic scales to relate the clocks 
on geodesics in bound systems, where space is not expand¬ 
ing, to the scale of comoving observers in the Hubble flow, 
where space is expanding. It is assumed without question 
that such clocks can be identified. 

In general relativity with the particular self-similar in¬ 
homogeneous geometry that arises from the evolution of 


the density perturbations of primordial inflation, the stan¬ 
dard solution of the fitting problem is not appropriate, as 
I argue in detail in a subsequent paper (Wiltshire 2005), 
henceforth Ref. I. In brief, while the underdense bubble, 
S, is on its largest scale a single perturbation away from 
the average bulk density, due to the scale-invariance of the 
inflationary spectrum, perturbations on the smaller and 
intermediate scales that give rise to the first structures 
nonetheless have the same statistical distribution as the 
bulk universe, with a mean density distributed about that 
of the bulk. The particular structure of inhomogeneity we 
see today, with its large voids, results in part from a “par¬ 
ticle horizon volume selection bias” in the initial density 
perturbations (see Ref. I). 

While the metric (2) is still valid for describing the 
clocks of an observer at an average spatial position, in 
an inhomogeneous underdense bubble such average spatial 
positions are in voids and do not coincide with stars and 
galaxies. On the other hand, the first bound systems of 
stars and star clusters which aggregate to typical galaxies, 
break away from the Hubble flow at early epochs when the 
average local density of matter in their past light cones is 
close to the average density of the bulk. Such systems 
retain fragments of the bulk hypersurface geometry. A 
clock in a system which has had an approximate station¬ 
ary Killing vector since the epoch of break away, thus mea¬ 
sures the bulk cosmic time parameter, r, which is “frozen 
in” from an earlier epoch in the evolution of the universe. 
Initial perturbations with the selection-biased underden¬ 
sity of the bubble, S, evolve to form voids. 

Given that structure forms in a bottom-up manner, a 
“temporal” Copernican principle operates within S in the 
sense that as average isotropic observers we assume we 
are located in a bound system formed from those density 
perturbations which broke away from the Hubble flow when 
the geometry within their past light cone was indistinguish¬ 
able from that of the bulk surfaces of homogeneity. Rather 
than occupying an average position on a spatial hypersur¬ 
face, we find ourselves in a bound system that formed from 
matter which broke from an almost homogeneous Hubble 
flow at an average epoch. At the present epoch, these 
same average galaxies, are located in clusters in bubble 
walls surrounding local voids, in a self-similar hierarchical 
structure, which we loosely call fractal. 

Since recent expansion of space occurs primarily in 
voids, the question of who does or does not measure an 
isotropic CMBR depends not only on local peculiar ve¬ 
locities but also on whether an observer’s line of sight 
to the surface of last scattering in any direction on the 
sky averages over the same number and volume distribu¬ 
tion of voids and bubble walls, once such structures form. 
For observers in average galaxies in bubble walls or ob¬ 
servers within small voids one would expect a roughly 
isotropic distribution of bubbles and voids within the self 
similar hierarchical matter distribution, and hence an al¬ 
most isotropic CMBR. Clock rates are determined by lo¬ 
cal geometry rather than by isotropy, or otherwise, of the 
CMBR at any location. 

To define the local cosmic time, r, in bound systems 
in average galaxies, as we have done, in reference to bulk 
hypersurfaces of matter homogeneity which stretch to re¬ 
gions beyond our present particle horizon may seem puz- 
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zling at first. However, prior to inflation, such regions 
were in causal contact with the observed universe. Thus 
the same basic processes of inflation that lead to isotropy 
of the CMBR also lead, via structure formation from scale- 
invariant perturbations, to a definition of inertial frames 
of observers in average galaxies, namely a new variant of 
Mach’s principle. (See Ref. I for further discussion.) 

Since general relativity is a local theory, the average ge¬ 
ometry in our first-order solution of the fitting problem is 
described by (2), (3), but insofar as measurements, includ¬ 
ing our own, are referred to the frame of average galaxies 
we must refer them to a coordinates related to those of 
( 2 ) by a non-trivial lapse function y(r) = so that ( 2 ) 
becomes 

ds 2 = 7 2 (r)ds 2 , (4) 

where 

ds 2 = —dr 2 + a 2 (r) + r 2 (cl 6> 2 + sin 2 6 d (j) 2 ) , (5) 

[ 1 + r 2 

and a = r )~ 1 a. Since quantities such as the luminosity dis¬ 
tance involve null geodesics, for many calculations we can 
use the conformally related geometry (5). It is not nec¬ 
essary to solve for the geodesics of observers, defined by 
local geometry. Only null geodesics probe the vast cosmo¬ 
logical scales over which the average geometry ( 2 ) is more 
relevant than other scales in the actual inhomogeneous 
geometry. Spacclike intervals are never directly measured. 
Thus by fiat, once we insist on writing the average cosmo¬ 
logical geometry in a “synchronous gauge” adapted to the 
proper time, r, of the local clocks of observers in average 
galaxies, then the geometry (5) also comes to define our 
rods when specifying areas, volumes and densities on the 
largest scales. The quantity j(t) leads to a gravitational 
time dilation between average galaxies in the gravitational 
wells of the bubble walls and the ideal comoving clocks in 
empty voids, as further discussed in Ref I. 

To analyse some quantities will require a specification of 
the entire inhomogeneous fractal geometry within S. This 
may be possible in the context of a LTB model. How¬ 
ever, given the hierarchical nature of the inhonrogeneity 
the LTB mass function would be much more complicated 
than those of the simple single- or few-void LTB mod¬ 
els that have typically been studied to date. Initially, we 
are only interested in properties defined by the average 
Hubble flow, although we might expect some modifica¬ 
tions of cosmological parameters due to the different range 
of scales and expansion rates among local voids (Tomita 
2001). While solving for intermediate scales, perhaps in 
terms of a LTB model, will ultimately be necessary, for 
the purpose of the first approximation in the fitting prob¬ 
lem, it suffices to construct a “spherical expansion model” 
of the underdense region S , in parallel to the well-known 
“spherical collapse model” (Kolb & Turner 1990). 

3. SPHERICAL EXPANSION MODEL 

Following the standard approach, we assume the initial 
density parameter is set sufficiently early that it is very 
close to unity: Hj = 1 — Si, 0 < Si -C 1. Furthermore, 
since the universe has flj ~ 1 initially we require the co¬ 
moving scales to match at that epoch. This means we set 
di ^ ai, 7 i ~ 1 , U ~ Ti and H t ~ 2/(3 n). 


There are various possible Hubble parameters to take 
into account: 

(i) the background Hubble parameter H = 2/[3(1 + ic)t] 
of the metric ( 1 ) which cannot be directly measured in S ; 

(ii) the Hubble parameter of comoving observers in (2), 
measurable only by unbound observers in average voids 

~, . 1 da 

H (t) = ; (6) 

a at 

(iii) the physical Hubble parameter we measure 

H(r) S if ; (7) 

(iv) an effective (unmeasured) parameter 

s If , (8) 

which is useful for the purposes of calculations. 

By our assumptions, while all these parameters coincide 
at early times, as the universe expands they may differ. 
As observers in average galaxies, using null geodesics we 
measure the cosmological geometry (5) referred to global 
cosmic time, r, not to the parameter t. We simply need 
to determine t(r) to determine cosmological quantities. 

Using Eqs. (2)-(5), the bulk universe scale factor, a = 
di (J^Hir ) n , where n = 2/[3(l + ie)] and Hi = 3nHi/2, the 
density parameter is defined as 

n = s,(5) 3 (M 2, AlS!AiA)iii. (#) 

\aJ \aij H 2 (cosh ?7 — l ) 3 

Comparing this with the standard parametric expression 
U(? 7 ) = 2(cosh?? — l)/sinh 2 77 derived from (2), (3), we 
obtain a parametric relation for r(rf), 

H IT = . (10) 

3(1 — Hi ) 3 / 2 sinhry 

The effective parameter ( 8 ) is found to be 
H H e ff p(l — H 0 ) 3 / 2 (H 0 + 2 )(cosh 77 + l ) 3 / 2 

H 0 (cosh 77 — l ) 3 / 2 (cosh 77 + 2 ) 

where fl Q = 2/(1 + cosh?? 0 ), and a subscript zero refers 
to the present epoch. The physically measured Hubble 
parameter (7) is 

H( x _ (cosh 2 ?; + 2 cosh 77 + 3)H ee (y) . . 

(cosh 77 + 1 ) (cosh 77 + 2 ) 

and the measured Hubble constant, H Q , is related to iL e ff 0 


/ a i\ 3 { a 

\d J \di 


Hr = 


H eS (v) = 


2 + SiJ 


Observe that at early times, 77 ~0, H e s~H and H~H, 
as expected but at late times, H | H. The lapse 

function is given by 

= dt = Hen 3 (cosh 77 + 1) 

1KI> dr h 2 (coshr 7 + 2 )‘ 1 ’ 

At the present epoch 7 0 = 3/(2 + H 0 ). 

We must also be careful to note that the locally mea¬ 
sured density parameter, differs from (9) by a volume fac¬ 
tor according to fi(r) = 7 3 (r)H(r), so that the present 
epoch measured matter density fraction is 

27H 0 

m (2 + H 0 ) 3 ’ 


(15) 
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with inverse = 6 f lm^ 2 sin 


,-i 


(Of) 


— 2 . 


The expansion age (10) is usefully rewritten as 

Q 0 (2 + ft 2 )(cosh T)~ l) 3 / 2 

r(r?) =-—---. 

H 0 ( 1 - n 0 )3/ 2 (n 0 + 2 ) 2 (cosh 77 + l ) 1 / 2 

The physically measured deceleration parameter q(r) = 
a/ 


(16) 


—H 2 a/a = —1 — H/H 2 , is given by 


q(v) = 


7 cosh 2 r] + 10 cosh 77 + I 
(cosh 2 77 + 2 cosh 77 + 3 ) 2 


(17) 


It is equal to the bulk deceleration parameter <1 = \ at 
early times 77 = 0 , but at late times as 77 —> 00 , q —> 0 , so 
that it is small but positive at the present epoch. This is 
just as expected from a model without dark energy, and 
contradicts the claims of cosmic acceleration by Kolb et 
al. (2005). However, we must recall that the supernovae 
measurements involve luminosity distances, and the inter¬ 
pretation of cosmic acceleration depends on the time pa¬ 
rameter assumed in taking derivatives. The present model 
may be recognised as mimicking a Milne universe at nearby 
redshifts, and this still provides a reasonable fit to present 
Snela data (Carter et al. 2005). 


4. OBSERVATIONAL TESTS 


To compare with observations Eqs. (12)-(17) need to be 
expressed in terms of the observed cosmological redshift, 
2 . Great care is needed at this point in identifying phys¬ 
ical quantities. It proves simplest to always use the line 
element (4), (5) remembering that r corresponds to clock 
time for observers in average galaxies. It follows that 

1 + z = % = |o2 = q-»o)( 2 +»o)( c ° s h^+i) _ (18) 
a a 7o H 0 (cosli ?7 + 2 )(cosh 77 — 1 ) 

The physical solution to eqn. (18) is 


cosh 77 


-1 

T 


(1 - n 0 )(2 + n p ) 
2fi 0 (* + l) 


0V[9f \z - 2H 2 + 16H 0 + 4] + (H 2 + 2) 2 

20p(;Z + 1) 

(19) 

We are now ready for our first cosmological tests. The 
luminosity distance is readily computed in the standard 
fashion and is found to be d L = (1 + 2 ) 7 0 _ 1 a 0 sinh (?7 — i] Q ), 

and since -v/ 1 — = 7 0 /(a 0 -ff e ff 0 ) it follows that 


As shown by Carter et al. (2005) the model gives a good 
fit to the supernova data, with or without non-baryonic 
cold dark matter in addition to baryonic matter. 

The expansion age is given by substituting (19) in (16). 
It gives a present day age of the universe of 


2(2+ ng) 

(2 + ^o ) 2 H 0 ' 


( 21 ) 


Using the best-fit value (Carter et al. 2005) for H Q = 
62.71-^7 km sec -1 Mpc -1 , the age of the universe is r Q = 
15.0 ± 0.4 Gyr for = 0.25 ± 0.5, or alternatively 
t q = 15.41q'| Gyr for a universe with only baryonic mat¬ 
ter, O m = 0.075 ± 0.055, using the recalibrated primordial 
nucleosynthesis bounds of Ref. I. 

Importantly, the expansion age is significantly larger 
at large 2 : with H Q = 62.7+1.7 km sec _1 Mpc^ 1 , for 
Q m = 0.25 ± 0.5 we find r = 4.3 ± 0.3 Gyr at 2 = 2, 
r = 1.451 q;) 2 Gyr at 2 = 6, and r = 0.31 ± 0.04 Gyr at 
2 = 20. For fl m = 0.075 ± 0.055 then r = 4.9^; 2 Gyr 
at 2 = 2, r = 1.94jjg;2i Gyr at 2 = 6, and r = 0.511 quo 
G yr at 2 = 20. This would buy precious time for struc¬ 
ture formation early on, and help to explain the redshift 
of the reionization epoch (Bennett et al. 2003), and obser¬ 
vations of “early” formation of structure (Cimatti et al. 
2004, Glazebrook et al. 2004). 


5. CONCLUSION 

As a model cosmology, the present Fractal Bubble Uni¬ 
verse is remarkable in that its gross features depend only 
on two already observed parameters, and the results ob¬ 
tained thus far are reasonably promising. If correct, then 
all observed quantities in cosmology must be systemati¬ 
cally reanalysed. The determination of some quantities 
will require further development of our understanding of 
the inhomogeneous matter distribution within S , including 
the largest scale local voids (Tornita 2001), in particular. 
However, the present model offers a clear framework for 
calculations, as is demonstrated in Ref. I, where further 
steps are taken in the reanalysis of the hot Big Bang and 
observational implications for the CMBR. 

It would be ironic that Einstein’s idea concerning the 
irrelevance of the cosmological constant, and his idea of 
the importance of Mach’s principle, may prove to both be 
right in understanding the universe. 




(1 + z)( 2 + n 2 ) 
fl 0 (2 + H q ) 


2 cosh 77 


(2 ^2q) 
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